Excited Heavy Quarkonium Production at the LHC through W^-Boson Decays 
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Sizable amount of heavy-quarkonium events can be produced through IF-boson decays at the 
LHC. Such channels will provide a suitable platform to study the heavy-quarkonium properties. 
The "improved trace technology" , which disposes the amplitude M at the amplitude-level, is helpful 
for deriving compact analytical results for complex processes. As an important new application, 
in addition to the production of the lower-level Fock states |(QQ')[1 , S']} an d I(QQ')[1-P])j we make 
a further study on the production of higher-excited \(QQ')) -quarkonium Fock states \(QQ')[2S]) , 
\{QQ')[3S]} and \ (QQ')[2P]} . Here \(QQ')) stands for the | (cc) }-charmonium, | (cfe)}-quarkonium and 
| (6fe))-bottomonium respectively. We show that sizable amount of events for those higher-excited 
states can also be produced at the LHC. Therefore, we need to take them into consideration for a 
sound estimation. If assuming all excited heavy-quarkonium states decay to the ground spin-singlet 
[l 1 So]-wave state with 100% efficiency via electromagnetic or hadronic interactions, we obtain the 
total decay width by adding all the mentioned Fock states together; i.e. by taking the bound-state 
parameters under the Buchmiiller-Tye potential model, we obtain ^w+-¥\(cc))+c3 ~ 591.3 KeV, 
IV+_>.|( c 5)>+!,s = 27.0 KeV, r w +->\(cl))+ce = 2 - 01 KeV . and r w+->\(bl))+cl = 93.3 eV. At the LHC 
with the luminosity C oc 10 34 cm _2 s _1 and the center-of-mass energy \/S = 14 TeV, this shows 



that 8.7 x 10 b 7] c and J/*, 4.3 x 10 5 B 
through W + decays. 



and £>*, 1.4 x 10 rjb and T events per year can be obtained 



PACS numbers: 12.38.Bx, 14.40.Pq, 14.70. Fm 



I. INTRODUCTION 

Within the framework of Non-Relativistic QCD 
(NRQCD) [J, a doubly heavy meson is considered as 
an expansion of Fock states. A systematic study on the 
production of lower |(Q<2'))-quarkonium Fock states at 
the 15-level {[(QQ'^l 1 S }) and \_(QQ')[1 S Si})) and the 
IP-level (KQQOtl'-Pi]) and \(QQ')[l 3 Pj])(J - 0,1,2)) 
via the W + semi-inclusive decays has been done in 
Ref.0. Here \(QQ')) 

-quarkonium stands for the |(cc))- 
charmonium, |(c6))-quarkonium and |(66))-bottomonium 
respectively. At the LHC, due to its high collision 
energy and high luminosity, sizable amount of heavy- 
quarkonium events can be produced through W + decays 
0, [I| ■ So these channels maybe an important supplement 
for other measurements at the LHC. 

In this paper, we will make a further discussion on 
the production of even higher |(Q(3'))-quarkonium Fock 
states via the W + semi-inclusive decays, which are at 
the nS-level (\(QQ')[nS})) and nP-level f](QQ')[nP])) 
with n > 2 respectively. As shown in Rcf.[2|, due to the 
color-suppression of the amplitude and the relative ve- 
locity suppression of the color-octet matrix element, the 
color-octet (QQ')" com P onei rt provides negligible contri- 
butions, so at the present, we will only discuss the color- 
singlet states' production channels. Moreover, because 
the W decay-widths to the |(&&))-bottomonium states 
are quite small, we will concentrate our attention on 
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the |(cc))-charmonium and the |(c6)) -quarkonium pro- 
duction. 

For the purpose, we will deal with the 1 — > 3 de- 
cay channel, W + \(QQ')[n]) + q + q' , which is a 
short notation for the following four semi-inclusive de- 
cay channels: W + — ► |(cc)[n]) +cs, — > \(cb)_[n]) +bs, 
W + -> \(cb)[n}) + cc and W + -> \(bb)[n}) + cb, respec- 
tively. Here n stands for the corresponding quantum 
number of the heavy KQQ'^-quat'konium Fock state, q 
and q' stand for the outgoing quark and anti-quark. In- 
tuitively, the heavy-quarkonium production could be un- 
derstood in terms of two distinct steps: one is the pro- 
duction of the heavy (QQ')-pair with certain quantum 
number n, together with the out-going quark q and anti- 
quark q'; and the second subsequent step is the evolution 
of such QQ' pair into a quarkonium bound state. The 
first step is pQCD calculable, where the quantum state 
[n] of the (QQ')-pa,ir can be achieved by a suitable pro- 
jector and the second step is characterized by a 
non-pcrturbativc matrix element which is proportional 
to the inclusive transition probability of the perturba- 
tive (<5Q')"Pai r with certain quantum number [n] into a 
bound state \(QQ')[n}) Q. 

To deal with the heavy-quarkonium production 
through the W + semi-inclusive decays, one needs to de- 
rive the pQCD calculable squared amplitude. The usual 
way is to deal with the squared amplitude |.M| 2 . Be- 
cause of the emergence of massive-fermion lines, the an- 
alytical expression for the squared amplitude becomes 
too complex and lengthy for more (massive) particles in 
the final states and for higher-level Fock states to be 
generated. For example, to derive the amplitudes for 



2 



the P-wave states, one also needs to get the derivative 
of the amplitudes over the relative momentum of the 
constitute quarks. It has been found that to do the 
numerical calculation using this conventional squared- 
amplitude technology becomes time-consuming for those 
complex processes, since the cross-terms of the matrix 
elements increase with the increment of Feynman dia- 
grams, \M\ 2 = J2ij ■M.i-Mj, where i and j stand for the 
number of Feynman diagrams of the process. 

One important way to solve this is to deal with the pro- 
cess directly at the amplitude level. For the purpose, the 
"improved trace technology" is suggested and developed 
in the literature 04HI- After generating proper phase- 
space points, one first calculate the numerical value for 
the amplitudes Aii, and then sum these values alge- 
braically and square it to get the squared amplitude, 
|A4| 2 = ISi-^i| 2 ; through such way, numerical sim- 
ulation efficiency can be greatly improved in compari- 
son to the conventional squared-amplitude technology. 
Moreover, under the approach, many simplifications can 
be done at the amplitude level due to the fermion-linc 
symmetries and the specific properties of each heavy- 
quarkonium Fock states, then, we can even written down 
the analytic expressions for the amplitude. 

The remaining parts of this paper are organized as fol- 
lows: In Sec. II, we give a short review of dealing the W + 
semi-inclusive decays by using the "improved trace tech- 
nology" . In Sec. Ill, we present the numerical results and 
discuss on the properties of the heavy-quarkonium pro- 
duction through W + decays. The final section is reserved 
for a summary. 



II. A SHORT REVIEW OF THE CALCULATION 
TECHNOLOGY 

Our present method for deriving the decay width of 
W+(k) -> \(QQ')[n})(q 3 ) + q(q 2 ) + q'(qi), especially the 
pQCD calculable part, is the same as that of Ref.[2|, for 
self-consistency, we will list the main points here. 

According to the NRQCD factorization formula (l2l |. 
the differential decay width of W + (k) \(QQ')[n])(q 3 ) + 
Q (92) + <f(<7i) can be factorized as 

dT = J2dt(W+ -> (QQ')[n] + qq")(0 H (n)) 7 (1) 
n 

where (O h (n)) describes the hadronization of a QQ' '-pair 
into the observable quark state H and is proportional 
to the transition probability of the perturbativc state 
{QQ')[n] into the bound state \{QQ')[n]) . The parame- 
ters k and qi are momenta of the corresponding out-going 
particles. 

The short-distance differential decay width 



df(W+ ->■ {QQ')[n] +qq' 



2k° 



J2\M\ 2 d$ 3 



(2) 



where means that we need to average over the spin 
states of the initial particles and to sum over the color 



and spin of all the final particles. In the W + rest frame, 
the three-particle phase space can be written as 
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The 1 — > 3 phase space with massive quark/ antiqark in 
the final state can be found in RefsJtJ, . With the help 
of the formulas listed in Refs.JH, [lO], one can not only 
derive the whole decay width but also obtain the corre- 
sponding differential decay widths that are helpful for ex- 
perimental studies, such as dT/dsi, dT/ds2, dT/dcos9i3 
and dr/dcos023, where si = (qi + qz) 2 , s 2 = (qi + q 2 ) 2 , 
#13 is the angle between qi and 03 in the W + rest frame, 
and 623 is the angle between q 2 and q 3 in the W + rest 
frame. Especially, the partial decay-widths over s\ and 
S2 can be expressed as: 



dT 



dsid,S2 



2567r 3 m^ / 



(4) 



where myy stands for the W-boson mass. The integration 
over si and so can be done with the help of the VEGAS 
program [l3| 0. 

Selection of appropriate angular momentum quantum 
number for the quarkonium state is done by suitable pro- 
jector and by performing proper polarization sum 0- 
|6(. The color-singlet non-perturbative matrix-element 
(Oi(n)) can be related to the Schrodinger wavefunctions 
at the origin |^(qq') (0) | or the first derivative of the wave- 
functions at the origin |VVqq')(0)I- For the color-singlet 
[nS]- and [nP\- wave states, we have 0, [l6| 



(O?(nS))^\^ im , )[nS]) (0)\ 2 , 
(O¥(nP))^^' im , )[nP]) (0)\ 2 . 



(5) 
(G) 



Here we have adopted the convention of Ref.Q for 
the non-perturbative matrix elements. Since the spin- 
splitting effects arc small, we will not distinguish the 
difference between the wavefunction parameters for the 
spin-singlet and spin-triplet states at the same n-th level. 



W+(k) 




(ceo [«])(«) 



W+(k) 
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Q'(m) 
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FIG. 1. Feynman diagrams for the process W + (k) — > 
\(cQ')[n])(q3) + Q' (q2)s(qi), where Q' stands for c or b quark 
accordingly, and |(cQ')[n]) stands for a heavy-quarkonium 
Fock state. 



1 The improved VEGAS version can be found in the programs 
BCVEGPY H and GENXICC QJ]. 
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FIG. 2. Feynman diagrams for the process W + (k) — > 
\(Qb)[n])(q3) + c(q2)Q(qi), where Q stands for the c or b quark 
accordingly, and |(Q6)[n]} stands for a heavy-quarkonium 
Fock state. 



The amplitude Ai for the mentioned channels can be 
calculated from two types of Feynman diagrams which 
are shown in Figs. (|1I2|) respectively. We present the two 
process W + — > \(cc)[n]) + cs and W + ~ > \(cb)[n]) + bs as 
W+(k) -> \{cQ')[n]){qi) + Q'(q 2 )s{qi), where Q' stands 
for c or b quark accordingly. We present the two processes 
W+_^ \(cb)[n])+cc and W+ -> \ (bb) [n])+cb as W+(k) 
\(Qb)[n])(q 3 ) + c(q2)Q{qi), where Q stands for the c or 
b quark accordingly. In those Feynman diagrams, the 
intermediate gluon should be hard enough to produce a 
cc pair or bb pair, so the amplitude is perturbative QCD 
calculable. 

For each production channel, the pQCD calculable am- 
plitude for a specific spin-state combination A4 SS > can be 
generally expressed as 



Cu, 



(7) 



where s and s' are spin indices, i and j are color in- 
dices for the outgoing quark and antiquark. The overall 

factor C = 2g £ v % KM Sh. Vckm stands for the Cabibbo- 

3 v6 J 

Kobayashi-Maskawa matrix clement, Vckm — V cs for 
W+ -> \(cc)[n]) +_cs and W+ -> \(cb)[n]) +_bs; Vckm = 
V cb for W+ -> \(cb)[n]) + cc and W+ -> \(bb)[n\) + cb. 

The expressions for A n can be read from Figs. (|1I2|) . 
and then we deal with the amplitude M by using the "im- 
proved trace technology" |7Hll1|. Under such approach, 
we first arrange the amplitude into four orthogonal sub- 
amplitudes M.± s ± s i according to the spins of the out- 
going quark q with spin s and antiquark q' with spin 
s'. After summing up the spin states of the outgoing 
quark/antiquark, the squared amplitude can be divided 
into four orthogonal parts, 



\M\ 2 



\M,\ 2 



\M 2 \ 2 



|X 3 | 2 + |^ 4 | 



(8) 



where the four amplitudes M.i (i = 1, • • • , 4) are defined 
as 

Mss'+M-s-s' Mss'-M-s-s' 

Mi = = , M 2 = — 



V2 

i a M s - s >-M- ss > 

= » M A = 



V2 

M s -s>+M- 



V2 



.(9) 



These four amplitudes M.i can be transformed into trace 
forms by properly dealing with the massive spinors with 



the help of an arbitrary light-like momentum kg and an 
arbitrary space- like momentum fci [k\ = —1), k and k\ 
satisfy the constraint ko ■ k\ = 0. 

More explicitly, we first introduce a massless spinor 
with negative helicity u_(fco) which satisfies the following 
projection 



U-(ko)u-(ko) = lu- fto, 



(10) 



where cj_ = (1 — "f 5 )/2. Then, we construct the positive 
helicity state as 



iU-(feo). 



It is easy to check that u + (fco)w+(fco) = u + where 
lu + = (1 + 7s)/2. Using these two massless spinors, the 
massive spinors for the fermion and antifermion can be 
written as follows: 



u s{q) = {4+ m)u^(k Q )/ y / 2k ■ q, 
u- s (q) = (4 + m)u + {k )/y/2k Q ■ q, 

v s(q) = (4- m)u-(ko) j\j2kQ ■ q, 
v- s (q) = {4- m)u + (k )/^2k ■ q, 



(12) 
(13) 
(14) 
(15) 



where the spin vector s M = ^ 

Using the above relations, the amplitudes M.i can be 
conveniently transformed into trace forms. And then we 
do the trace of the Dirac 7-matrix strings which will re- 
sult in explicit series over some limited and independent 
Lorentz-structures. The final results are independent of 
the choice of fco and k\ , one can choose them to be those 
that can maximumly simply the analytical expressions 
of the amplitude. The compact and complex expressions 
for the amplitudes of the S'-wave and P-wave cases can 
be found in the appendix of Ref. @ • 

These results provide the foundation for estimating the 
higher Fock-states' properties. 

Even though it is hard to get the analytic results for the 
conventional squared-amplitude technology, as a cross- 
check, we also adopt it to calculate the VF + -decay widths 
numerically. Under the same parameter values, we have 
found agreement of these two approaches. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Input parameters 

When doing the numerical calculation, the iirput pa- 
rameters are chosen as the following values [13, 13 : 
m w = 80.399 GeV, T w + = 2.085 GeV, m B = 0.105 
GeV, |U CS | = 1.023 ± 0.036 and |U c6 | = 0.0406 ± 
0.0013. Leading-order a s running is adopted. An op- 
timal process to set the renormalization scale [i r has re- 
cently been suggested in the literature, i.e. the Prin- 
ciple of Maximum Conformality (PMC) [l^, which is 
rcnormalization-scheme independent and can provide al- 
most renormalization-scale independent estimation even 
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n = 1 


n = 2 


n = 3 


|-R|(cc)[nS]>(0) 2 
\R\(cc)\nP]) (0)1 


0.810 GeV 3 
0.075 GeV 5 


0.529 GeV 3 
0.102 GeV 5 


0.455 GeV 3 


|-R|(c5)[nS]>(°) P 
l W |( C b)[nP]}( U ) 


1.642 GeV 3 
0.201 GeV 


0.983 GeV 3 
0.264 GeV 


0.817 GeV 3 


-R|(bb)[nS]>(0) P 
|- R |(i>6)[nP])( )l 2 


6.477 GeV 3 
1.417 GeV 5 


3.234 GeV 3 
1.653 GeV 5 


2.474 GeV 3 


m c 


1.48 GeV 


1.70 GeV 


1.90 GeV 


nib 


4.88 GeV 


5.00 GeV 


5.10 GeV 



TABLE I. Bound-state parameters adopted in the calculation, 
which are derived under the Buchmuller-Tye potential [20l |. 



at the fixed-order. However to set the leading-order 
PMC scale, one needs to know the next-to-leading order 
{/?i}-terms which will be absorbed into the a s -running 
to form a commensurate leading-order PMC scale. For 
the present leading-order calculation, since T oc a 2 s {i± r ), 
the renormalization scale-uncertainty can be easily fig- 
ured out when we know the physical scale well. For clar- 
ity, we set [i r to the conventional choice of 2m c , and 
a s (2m c ) = 0.26. 

As for the wavefunction at the origin and the first 
derivative of the wavefunction at the origin, we adopt 
the values derived in Ref. [2(| under the Buchmiiller-Tye- 
potential (BT-potential) model as their central values, 
since it is noted that the BT-potential has the correct 
two-loop short-distance behavior in pQCD [2l|. The re- 
sults for three other potential models, i.e. the Power-Law 
model 12211 . the Logarithmic model [23| and the Cornell 
model [24j , will be adopted as an error analysis. Further- 
more, similar to our previous treatment @], we adopt the 
same constitute quark masses for the same n-th level Fock 
states. To ensure the gauge invariance of the hard am- 
plitude, we set the | (QQ'^-quarkonium mass M to be 
m Q + m Q'- 

We present the quarkonium bound-state parameters 
under the BT-potential model in Table HI where the 
constitute charm- and bottom- quark masses, the ra- 
dial wavefunction at the origin and the first derivative 
of the radial wavefunctions at the origin for \(QQ') [n])- 
quarkonium are presented. Here RuQQ>\i n g-i\(0) and 



1(QQ')[nS]>< 

l 1(QQ')[™- p ]) ^ are re ^ a ^ e< ^ t° the wavefunction at the ori- 
gin and the first derivative of the wavefunction at the 
origin through the following equations: 



\(QQ')[nS]) 
\(QQ')[nP]) 



(0)1 

(0)| = y^\R\ 



VV4^|-R|(QQ')[nS])(0)|, 
|(QQ')[nP]>(°)|- 



(16) 



B. Heavy-quarkonium production via W + decays 

The decay widths for the aforementioned quarko- 
nium states through the production channel, W + — ¥ 
|(Q<2')N) + 1Q' > are listed in Tables HJ and Hill 





n=l 


n=2 


n=3 


W + -» Kcc^So]) +cs 


132.0 


56.36 


22.85 


W + ->■ \(cc)[n a Si]) + cs 


136.4 


58.30 


33.38 


W + \(cb)[n 1 So\) +bs 


6.39 


3.53 


2.74 


W + -> \(cb)[n 3 Si]) +bs 


5.49 


3.07 


2.41 


W+ -> \(cb)[n 1 S ]} +cc 


0.411 


0.160 


0.094 


W+ -> |(c*)[ra 3 ii?i]) +cc 


0.593 


0.224 


0.128 



TABLE II. Decay widths (in unit: KeV) for the production 
of |(QQ')[nS])-quarkonium through W + decays. 





n=l 


n=2 


W+ -> 


(cc)[n 1 Pi]) + cs 


22.95 


12.98 


W+ -> 


(cc)[n 3 P ]) + cs 


28.33 


17.21 


W+ -»■ 


(cc)[71 3 Pi]} + cs 


28.31 


19.07 


W+ -»■ 


(cc)[n 3 P 2 ]) + cs 


14.76 


8.414 


W+ -»■ 


(c6)[n 1 Pi]) +bs 


0.270 


0.360 


W+ -> 


(cb)[n 3 P ]} +bs 


0.733 


0.693 


W+ -»■ 


(cb)[n 3 Pi]) +bs 


0.514 


0.724 


w + -> 


(cb)[n 3 P 2 ]) +bs 


0.034 


0.039 


W+ -> 


(cb^r^Pi]) + cc 


0.105 


0.066 


W+ -> 


(cb)[ji a P ]) + cc 


0.026 


0.019 


W+ -> 


(c6)[n 3 Pi]) + cc 


0.054 


0.036 


W+ -»■ 


(cb)[n 3 P 2 ]) + cc 


0.060 


0.037 



TABLE III. Decay widths (in unit: KeV) for the production 
of |(QQ / )[nP])-quarkonium through W decays. 



From Tables |H] and IIII1 it is found that in addition to 
the ground 15- level states, the higher (<5<5')-c[uarkoniuni 
states can also provide sizable contributions to the total 
decay width; i.e., 

• For the charmonium production channel, W + — > 
\(cc)[n])+cs, the decay widths for [n] = 25, 35, IP 
and 2P-wave states are about 43%, 21%, 35% and 
21% of that of the 15-level wave state, respectively 

• For the (cfr)-quarkonium production channel, 
W + \(cb)[n])+bs, the decay widths for [n] = 25, 
35, IP and 2P-wave states are about 55%, 43%, 
13% and 15% of that of the 15-level, respectively. 

• For the (c6)-quarkonium production channel, 
W + \(cb)[n]) +cc, the decay widths for [n] = 25, 
35, IP and 2P-wave states are about 38%, 22%, 
24% and 16% of that of the 15-level, respectively. 

Here, for convenience, we have used [nS] to represent the 
summed decay width of [n 1 5o] and [n 3 5i] at the same 
n-th level; [nP] to represent the summed decay width of 
[n x Pi] and [n 3 Pj] at the same n-th level. 

To show the relative importance among different Fock- 
states more clearly, we present the differential distribu- 
tions dT/dsi, dT/ds2, dT/dcos9i3 and dT/d cos #12 for 
the mentioned channels in Figs. (|3l4l5p and Figs. (|6l7l8|) . 



s,(GeV 2 ) 



s 2 (GsV z ) 



FIG. 3. Differential decay widths dT/dsi and dT/d,S2 for 
W + — > \{cc)[n\) + cs, where the diamond, the dash-dot, 
the dotted, the solid and the dashed lines are for |(cc)[lS]}, 
\(cc)[2S\), \{cc)[3S]), |(ec)[LP]) and \(ce)[2P]}, respectively. 



FIG. 7. Differential decay widths dr/Gfcos(9i3 and 
dr/d cos #23 for W + — > \{cb)[n\) + bs, where the diamond, 
the dash-dot, the dotted, the solid and the dashed lines are 
for \(cb)[lS]}, \(cb)[2S]), \(cb)[3S]), \(S)[1P]) and \(cb)[2P]), 
respectively. 




s (GeV a ) s (GeV 2 ) 



FIG. 4. Differential decay widths dT/dsi and dT/ds2 for 
W + — > \(cb)[n]) + bs, where the diamond, the dash-dot, 
the dotted, the solid and the dashed lines are for |(c5)[f5']), 
\(cb)[2S}), \(cb)[3S}), \(cb)[lP}) and \{cb)[2P}), respectively. 




COS9 13 COS<) 23 



FIG. 8. Differential decay widths cir/dcos6i3 and 
dr/dcos#23 for W + — > \{cb)[n]) + cc, where the diamond, 
the dash-dot, the dotted, the solid and the dashed lines are 
for |(c6)[lS]>, \(cb)[2S]}, \(cb)[3S\), \(cb)[lP]} and \(d>)[2P]), 
respectively. 





FIG. 5. Differential decay widths dT/dsi and dT/ds2 for 
W + — y \(cb)[n]} + cc, where the diamond, the dash-dot, 
the dotted, the solid and the dashed lines are for |(c&)[15']), 
\(cb)[2S]}, \(cb)[3S}}, \(cb)[lP]) and \ (cb)[2P}}, respectively. 




FIG. 6. Differential decay widths dr/a!cos#i3 and 
dr/dcos6 l 23 for W + — y \(cc)[n]} + cs, where the diamond, 
the dash-dot, the dotted, the solid and the dashed lines are 
for \(cc)[lS}}, \(cc)[2S\), \(cc){3S]}, \(cc)[lP]) and |(cc)[2P]>, 
respectively. 



Moreover, taking the channel W + 
example, we define a ratio 



|(cc)[n]> 



cs as an 



Ri[n] = 



dT/d Si (\(cc)[n})) 
dT/dsiiKcc^lS}))' 



(17) 



where i = 1,2. The curves are presented in Fig. ([9]). 
These figures show explicitly that higher Fock-states 
\(QQ')[2S}), \(QQ')[3S]), \(QQ')[1P}) and \{QQ')[2P}) 
can provide sizable contributions in comparison to the 
lower Fock-state KQQ')!!^]) in almost the whole kine- 
matical region. 

If assuming all the higher excited heavy-quarkonium 
states decay to the ground spin-singlet S'-wave state 
((QQ^Wl 1 So})) with 100% efficiency via electromagnetic 
or hadronic interactions, we obtain the total decay width 



r 



W+-H(cc))+cs 
W+^\(cb))+bs 



591.3 KeV, 
27.0 KeV, 



Fw+->\{cb))+cc — 2 -01 KeV, 



W+^\(bb))+cb 



93.3 eV. 



(18) 
(19) 
(20) 
(21) 



It shows that the total decay width of the bottomonium 
case is small as previously mentioned. 



At the LHC, with the luminosity C 



p-p 



cm 



and the center-of-mass energy y/s = 14 TeV, about 
3.07 x 10 10 of W + events per year can be produced 
[25|. Then, we can estimate the event numbers of 
|(QQ'))-quarkonium production through W + decays; i.e. 




S^GeV 2 ) s 2 (GeV 2 ) 



FIG. 9. The ratios Ri[n] and 7?2[n] versus si and S2 for the 
channel W + — > \(cc)[n]} + cs. Here the dotted, the dash-dot, 
the solid and the dashed lines are for |(cc)[2S]), Kcc)^.!?]), 
|(cc)[lP]) and |(cc)[2P]), respectively. 





B.T.[211 


P.L.[22] 


Log. [23] 


Cor.[24j 


[n] = [l 1 ^] 


132.0 


162.8 


132.8 


236.9 


[n] = [l 3 Si] 


136.4 


168.3 


137.3 


244.9 


[n] = [2 1 S'o] 


56.36 


59.55 


44.53 


98.76 


[n] = [2 3 Si] 


58.30 


61.60 


46.06 


102.2 


[n] = frSo] 


22.85 


20.59 


14.36 


39.72 


[n] = [3 3 ,Si] 


33.38 


30.08 


20.98 


58.04 


[n] = [l^i] 


22.95 


38.25 


23.87 


40.09 


M = [l 3 Po] 


28.33 


47.22 


29.47 


49.49 


[n] = [l 3 Pi] 


28.31 


47.18 


29.44 


49.45 


[n] = [1 3 P 2 ] 


14.76 


24.59 


15.35 


25.77 


[n] = [2^1] 


12.98 


16.67 


9.673 


23.67 


[n] = [2 3 P ] 


17.21 


22.10 


12.82 


31.38 


[n] = [2 3 Pi] 


19.07 


24.50 


14.21 


34.78 


M = [2 3p 2] 


8.414 


10.81 


6.269 


15.34 


sum 


591.3 


734.2 


537.1 


1050. 



TABLE IV. Decay widths (in unit: KeV) for (cc)-charmonium 
production channel W + — > \(cc)[n]} + cs, where bound-state 
parameters from four potential models are adopted. 



3.95 x 10 6 \(cc)[lS]), 1.69 x 10 6 \(cc)[2S]}, 8.29 x 10 5 
\(cc)[3S]), 1.39 x 10 6 |(cc)[lP]), 8.49 x 10 5 \(cc)[2P]) 
events per year can be produced; 1.90 x 10 5 |(c6)[15]), 
1.03 x 10 5 \(cb)[2S}), 7.91 x 10 4 \(cb)[3S}), 2.64 x 10 4 
\(cb)[lP]), 2.91 x 10 4 \(cb)[2P]) events per year can be 
produced. 



C. Decay widths under four potential models 

Next, we discuss the uncertainty caused by the 
bound-state parameters. These parameters are main 
uncertainty sources for estimating the heavy \(QQ'))- 
quarkonium production. We take the parameters derived 
under four potential models, i.e. the Buchmuller-Tye 
model [2l| , the Power-Law model 1221 , the Logarithmic 
model [231 an d the Cornell model [24| , to do our discus- 
sion. We take the symbols B.T., P.L., Log. and Cor. 
as short notations for the Buchmullcr-Tyc model, the 
Power-Law model, the Logarithmic model and the Cor- 
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B.T.[211 


P.L.[22] 


Log. [23] 


Cor.[24j 


[n] = [l^o] 
[n] = [l 3 Si] 


6.39 
5.49 


6.66 
5.72 


5.87 
5.05 


12.4 
10.7 


[n] = pSo] 
[n] = [2 3 Si] 


3.53 
3.07 


3.41 
2.97 


2.76 
2.40 


6.33 
5.51 


[n] = frSo] 
[n] = [3 3 Si] 


2.74 
2.41 


2.28 
2.01 


1.89 
1.66 


4.84 
4.26 


[n] = [l'Pi] 

N = [i 3 Po] 

[n] = [l 3 Pi] 
[n] = [1 3 P 2 ] 


0.270 
0.733 
0.514 
0.0336 


0.439 
1.19 
0.836 
0.0547 


0.321 
0.871 
0.611 
0.0400 


0.459 
1.25 
0.875 
0.0572 


N = l 2lp i] 
[n] = [2 3 P ] 
[n] = [2 3 Pi] 
M = [2 3 -P2] 


0.360 
0.693 
0.724 
0.0392 


0.480 
0.924 
0.965 
0.0522 


0.326 
0.627 
0.655 
0.0355 


0.629 
1.21 
1.26 

0.0684 


sum 


27.0 


28.0 


23.1 


49.8 



TABLE V. Decay widths (in unit: KeV) for (cfe)-quarkonium 
production through W + — > \ (cb)[n]} + bs, where bound-state 
parameters from four potential models are adopted. 





B.T.[211 


P.L.[22J 


Log.[23j 


Cor.[24j 


[n] = [l^o] 
[n] = [l 3 Si] 


0.411 
0.593 


0.428 
0.618 


0.378 
0.545 


0.798 
1.15 


[n] = [2^0] 
[n] = [2 3 Si] 


0.160 
0.224 


0.155 
0.216 


0.126 
0.175 


0.288 
0.402 


[n] = frSo] 
[n] = [3 3 Si] 


0.0942 
0.128 


0.0784 
0.107 


0.0649 
0.0883 


0.166 
0.227 


[n] = [l'Pi] 
N = [l 3 -Po] 
[n] = [l 3 Pi] 
[n] = [1 3 P 2 ] 


0.105 
0.0255 
0.0536 
0.0598 


0.170 
0.0415 
0.0872 
0.0972 


0.124 
0.0303 
0.0637 
0.0711 


0.178 
0.0434 
0.0912 
0.102 


[n] = [2 x Pi] 
[n] = [2 3 P ] 
[n] = [2 3 Pi] 
H = [2 3 -P2] 


0.0659 
0.0185 
0.0364 
0.0374 


0.0878 
0.0247 
0.0486 
0.0498 


0.0596 
0.0168 
0.0330 
0.0338 


0.115 
0.0324 
0.0636 
0.0653 


sum 


2.01 


2.21 


1.81 


3.72 



TABLE VI. Decay widths (in unit: KeV) for (cfe)-quarkonium 
production through W + — > \(cb)[n]) + cc, where bound-state 
parameters from four potential models are adopted. 



nell model, respectively The constitute quark masses 
and their corresponding radial wavefunctions at the ori- 
gin and the first derivative of the radial wavefunctions 
at the origin for the \(QQ'))-quar]iomxim states can be 
found in Tables I, II and III of Ref.[lo!]. To short the 
paper, we do not repeat them here. 

The decay width for | (QQ'))-qi\ai konium production 
under four potential models are presented in Tables 
IIVI IVl and IVII The Cornell model, which is a naive 
Coulomb-plus-linear potential, always gives the largest 
values among the four models. While, the decay widths 
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for the other three models are consistent with each other: 
taking the BT-model decay width as the center value, for 
the channel W + — > \(cc)[n]) + cs, we obtain the uncer- 
tainty where the upper value is for the Power- 
Law model and the lower value is for the Logarithmic 
model; for the channel W + — > \(cb)[n])+bs, we obtain the 

uncertainty (^5%) > ^ or ^ e cnann °l ^ + — ► l( c ^)[ n ]) + C A 
we obtain the uncertainty ±10%. 

IV. SUMMARY 

In addition to our previous studies on the lower Fock 
states' production presented in Ref. 0, we have car- 
ried out an investigation on the higher-excited \(QQ))- 
quarkonium states' production through lU + -boson semi- 
inclusive decays. 

For the | (cc) )-charmonium production channel W + — > 
\(cc)[n]) + cs, the decay width for [n] = 2S, 35, IP and 
2P states are about 43%, 21%, 35% and 21% of that of 
the 15-level; For the |(c6)) -quarkonium production chan- 
nel W + -> \(cb)[n}) + bs. the decay width for [n] = 25, 
35, IP and 2P states are about 55%, 43%, 13% and 15% 
of that of the 15-level; For the |(c6))-quarkonium pro- 



duction channel W + — > \(cb)[n]) + cc, the decay width 
for [n] = 25, 35, IP and 2P states are about 38%, 22%, 
24% and 16% of that of the 15-level. Then these higher- 
excited n5-level and nP-level states with n > 2 can also 
provide sizable contributions to the heavy quarkonium 
production. Therefore, we need to take these higher ex- 
cited states into consideration for a sound estimation. 

At the LHC with the luminosity C cx 10 34 cm _2 s _1 and 
the center-of-mass energy \fS = 14 TeV, with the help of 
Eqs. (|18|19|20|21J) . one can obtain 8.7 x 10 6 r] c and J/*, 
4.3 x 10 5 B c and B* and 1.4 x 10 3 rfc and T events per 
year. Uncertainties caused by the bound-state parame- 
ters arc presented, four potential models arc adopted for 
the purpose. It is found that the decay widths for the 
BT model, the Power-Law model and the Logarithmic 
are consistent with each other, whereas the naive Cornell 
model is not. 
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